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We study the decay B — > J/ijjK tt^-k using 117 million BB events collected at the y(4S') 
resonance with the BaBar detector at the PEP-II e^e~ asymmetric-energy storage ring. We mea- 
sure the branching fractions B {B~ -> J/iIjK'tv'^tt') = (116 ± 7{stat.) ± 9{syst.)) x 10"^ and 
B {B- -^ X{3S72)K-)x B (X(3872) ^ J/V'Tr+Tr") = (1.28 ±0.41) x lO"'"' and find the mass of 
the X(3872) to be 3873.4 ± 1.4MeV/c^ We search for the he narrow state in the decay B~ -> h^K^ , 
he -* J/ipn+TT- and for the decay B" -^ J/ijjD^-k-, with D° -^ K--K+. We set the 90% C.L. limits 
B {B' -> hcK') X B (hc^ J/^Tx+Tx-) < 3.4 x lO^** and B {B~ ^ J/i;D°n-) < 5.2 x 10"^ 

PACS numbers: 13.25.Hw, 12.15.Hh, ll.30.Er 



The study of B decays to final states containing 
charmonium and strange mesons is especially suited to 
the search for new charmonium states and for intrinsic 
charm. In particular, the decay B~ -^ J/tpK~TT^TT~ [IJ] 
can occur via the production of charmonium states de- 
caying into J/iJjTT^TT^ or possibly via B~ —> J/jpD'^i:~ , 
with D^ ^ K- TT+. Recently the Belle 3 and CDF 3] 
collaborations have observed a new state, the X(3872), 
decaying into J/^t:^-k^ . This state is a charmonium 
candidate, the l^Dj or l^Dg j||, with J^^ = 2 — or 
J^'-^ — 3 , or even possibly a molecule of charmed D 
and D* mesons |5(. In this Letter, using llTmillion T{AS) 
decays into BB pairs, we confirm the observation of the 
X(3872) and search for the unconfirmed charmonium IPi 
state /ic(3526) [6j. In addition, we study final states in- 
volving D mesons to test models developed to explain 
the excess of low momentum J/if) mesons in inclusive B 
decays [ij. 

The data were collected at the PEP-II asymmetric- 
energy e^e~ B-factory with the BABAR detector, which 
is fully described elsewhere |2|. The detector includes 
a silicon vertex tracker and a drift chamber in a 1.5- 
T solenoidal magnetic field, which detect charged par- 
ticles and measure their momentum and energy loss. 
Photons, electrons, and neutral hadrons are detected in 
a CsI(Tl)-crystal electromagnetic calorimeter. A ring- 
imaging Cherenkov detector is used for particle identifi- 
cation. Penetrating muons and neutral hadrons are iden- 
tified by resistive-plate chambers in the steel of the flux 
return. We use a Monte Carlo simulation of the BABAR 
detector based on GEANT4 [iflj to validate the analysis 
procedure and to estimate efficiency corrections. 

The event reconstruction and selection follow closely 
those described in an earlier paper [ll|. The present 
analysis has been optimized to maximize the sensi- 
tivity to B^ -^ J/tpK~TT~^Tr~ decays. We reconstruct 
J/ip -^ e+e~ candidates from pairs of tracks selected with 
criteria that are 98% (7%) efficient for electrons (pions). 
To account for energy losses, we combine the electron 
pairs with bremsstrahlung-photon candidates and use an 
asymmetric mass window, 2.95 < m.ee(7) < 3.14GeV/c^. 
We reconstruct J/ip — > ^-^ It-^ candidates from pairs of 
tracks selected with criteria that are 77% (8%) efficient 
for muons (pions), satisfying 3.06 < m^^ < 3.14GeV/c^. 



The nominal J/tp mass |l2| is imposed as a constraint 
on J/f/; candidates, thereby improving the resolution on 
the B mass and on any charmonium states in its de- 
cay. Kaons are identified using criteria that have an ef- 
ficiency of 97%, with a 15% pion-misidentification rate. 
B-meson candidates are formed by combining a J/tp can- 
didate with a kaon candidate and two additional op- 
positely charged tracks. To suppress further the back- 
ground from light-quark production, which is character- 
ized by back-to-back jets, the angle 9t between the thrust 
axes of the reconstructed B candidate and the rest of the 
event in the center-of-mass system is required to satisfy 
|cos6't| < 0.8(0.9) for J/tp -> e+e~ {J/tJj -^ fi+fi^) can- 
didates. 

Signal and combinatorial background are discrimi- 
nated using two kinematic variables: the beam-energy- 
substituted mass, toes = \/(Vs/2)^ —p|j^, and the dif- 
ference of the B candidate's measured energy from the 
beam energy, AE = Eg — {y^/2). Here Eg {p*g) is the 
energy (momentum) of the B candidate in the center- 
of-mass frame and \/s is the total center-of-mass energy. 
The signal region is defined to be |AB| < 3ct, where the 
resolution ct, determined with data, is 12MeV. A binned 
likelihood fit to the toes distribution (Fig.QIa)) is used 
to separate the signal, taken as a Gaussian distribution 
with a fitted width of about 2.5MeV/c^, plus a small tail 
to account for energy losses |l4| . from the combinatorial 
background distributed as an ARGUS threshold func- 
tion [ig]. We have checked with Monte Carlo simulation 
that there is no significant background from B decays 
that has the same jties distribution as the signal. 

To reduce systematic uncertainties, we measure 
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where Nevents — 2540 ± 72 is the number of 
B^ -^ J/ipK^ TT^TT~ signal events extracted from the fit 
to the Toes distribution. The number of 'ip(2S) events, 
-^i/j(2S) — 556 ± 30, is obtained by fitting mj^^^T^ dis- 
tribution, after subratcting combinatorial background, 
with two Gaussian distributions representing the ^{2S) 
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FIG. 1: Distribution of rriES for (a) B~ —> J/ip K~-k'^-k~ 
candidates, and (b) events in the X(3872) region, 3862 < 
mj/ipivTv < 3882 MeV/c . The solid curves represent the 
binned likelihood fits described in the text; the combinato- 
rial components are indicated by the dashed curves. 



FIG. 2: Distribution of mj/.^.^.^ (a) in the entire range, (b) 
in the he region, (c) at the ^/>(2S), and (d) in the region of 
the X(3872) with the projection of the unbinned likelihood 



fit superimposed, 
applied. 



The requirement rriES > 5.27 GeV/c is 



signal and a flat distribution representing the remaining 
background. (Fig. \^c) shows the corresponding unsub- 
tracted distribution). Throughout this Letter the distri- 
butions after combinatorial-background subtraction are 
obtained by fitting the tjies distribution of the events 
within each bin of the variable of interest (?7i j/^tttt i^i this 
case). The binned x^ fit gives a resolution on nij/^T^^ of 
3.1 ± 0.2 MeV/c^ for the core Gaussian containing 70% of 
the events and 12 ± 3MeV/c^ for the broader Gaussian. 
The total B~ -^ J/tpK~TT~^TT~and the 'ip{2S) selection ef- 



ficiencies, e and e 



V(2S) 



are extracted from Monte Carlo 



simulation: we obtain e^(2S)/e = 1.17 ± 0.03. We use 
Bi^i2S) -^ J/ip-K+n-) = (31.8 ±1.0)% H^- 

We estimate the systematic error due to the choice 
of the signal m^s shape function by replacing it with a 
simple Gaussian. We estimate the uncertainty on the fit 
to the rrij/^^^ distribution by using the signal resolution 
function as measured on Monte Carlo and by varying the 
background shape. Including all these errors, we measure 
R = 1.70±0.10(stat.)±0.09(s2/st.) which, combined with 
B{B- -^ V(2S') K-) = (6.8 ± 0.4) x lO"'* [13, yields 

B{B- ^J/ijK-TT+TT-) = (2) 

(116 ± 7{stat.)±9{syst.))xl0-^. 

To investigate the possible presence of narrow char- 
monium states decaying to J/ip '!t~ n'^ , we have studied 

(Fig. |2Ia)). We observe an 



excess in the region of the X(3872) (Fig. I^Jd)), but do 
not find any excess in the he region (Fig. [2Ib)). The 
mass of the X(3872) state is extracted from an un- 
binned maximum likelihood fit to the two-dimensional 
distribution in mss Siiid rrij/^^^ . The probability den- 
sity function (PDF) is taken to be the sum of four 
terms. The first three describe B^ -^ J/ipK^-K^TT^ de- 
cays that peak when mss is the mass of the S-meson. 
The PDF of these three terms contains a Gaussian func- 



tion in rriEs times a function of m 



J/^TT7 



that describes: 



1) non-resonant events, distributed as a first order poly- 
nomial; 2) ip{2S) candidates, distributed as a double- 
Gaussian resolution function around a mean value that 
is allowed to fioat; and 3) X(3872) candidates, with the 
same resolution function as the tpi^S) but with a mass 
that floats relative to the 'ip{2S) mass. The measure- 
ment of mass difference allows us to neglect systematic 
errors on the absolute mass scale. The fourth term of 
the PDF describes the combinatorial background, dis- 
tributed as an ARGUS threshold function in toes and 
as a first order polynomial in rnj/^T^j^. From the i^{2S) 
mass value, to^(2S) = 3685.96 ± 0.09 MeV/c^ [13, we find 
77ix(3872) = 3873. 4± 1.4MeV/c^, consistent with the pre- 
vious measurements by Belle 2j and CDF 3] . 

The measurement of the branching fraction 



the distribution in m 



J/l/fTl 



B{B- -^ X{3872)K-) x B{X{3872) - 
performed with a counting technique. 



J/^TT'^n ) is 
We select events 



in a ±10MeV/c^ window around mj/ ,^^^=3872 MeV/c^, 
and find the number of events with tties > 5.27 GeV/c^ to 
be N^ata = 63. We estimate the number of these events 
due to combinatorial background {Ncomb — 22.0 ± 4.3) 
from a fit to the toes distribution (Fig. db)). 
The number of events with the same final state 
B- -^ J/il^K-TT+TT-, but not belonging to the X(3872) 
signal, is estimated to be Npeak = 10.5 ± 3.2 from a 
fit to the TO.ES distribution in the symmetric sideband 
15 < \mj/^^^ - 3872| < 45MeV/c^ The resulting 
number of signal events is 30.5 which agrees within 
the errors with the number of signal events, 25.4 ± 8.7, 
obtained from the fit to the X(3872) in Figure [^d). The 
branching fractions are determined using a frequentist 
confidence level jl3| . This technique treats properly 
the small number of events and includes the system- 
atic errors directly in the computation of confidence 
intervals or limits. The confidence level, a, a function 
of B{B- -. X{3872)K-) x B{X{3872) -^ J/^tt+tt-) is 
computed as the fraction of times that a random number 
generated according to a Poisson distribution with a 
mean value of 

A< = Nbkg + N^(^2S)f-w (3) 

B{B- -^ X{3872)K-)B{X{3872) -^ J/i/'tt+tt-) 
B{B- -^ i^{2S)K-)B{tlj{2S) -^ J/V'tt+tt-) 

exceeds the observed data. For a given value of 
BiB- -^ XiS872)K-) x e(X(3872) ^ J/V'tt+tt-) the 
variables Nbkg, N^{2S), B{B^ -^ V(2S')i4:"), and 
B{'tp{2S) -^ J/-07r+7r~) are randomly generated to de- 
termine a value of /i, which is then used in a Poisson 
distribution to generate a new value of the number 
of detected events. The generation is repeated many 
times and the fraction of times the random num- 
ber exceeds Ndata = 63 yields the value of a. The 



variables Ni 



bkg: 



N, 



^(2S): 



B{B' 



^i2S)K-) 



and 



B{tfj{2S) -^ J/V'TT+TT ), are generated according to 
Gaussian distributions. The mean of N^(^2S) is 556 and 
a = 30. The mean of Nbkg is Ncomb +Npeak = 32.5 
and a — 5.9, which includes a systematic error on Npeak 
calculated by varying the boundaries of the sideband. 
We use published values |12i] for the remaining branch- 
ing fractions and their errors, assumed to be Gaussian. 
Finally, e^ = (92 ± 1)% is the fraction of events that 
fall in the mj/^^j^ window, from applying the same 
mass window cut to the ip{2S) and assuming the same 
efficiency. From the values of B{B- -^ X{3872)K~) at 
which a = 16% and 84% we measure 



B{B- -^ X{3872)K- 



X B{X{3872) -^ J/tPn+TT-) = 
(1.28 ± 0.41) X 10"^ (4) 



The probability that the observed events are a 
background fluctuation in the considered mass 
window is 5.4 x 10"'*, corresponding to 3.5 Gaus- 
sian standard deviations. As a check, we per- 
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FIG. 3: Distribution of m^+^- (a) at the X(3872) and (b) 
at the tf}{2S), after subtraction of combinatorial and peaking 
background. 



formed the same measurement on the J/i/j — *■ e+e 
and J/ip — * H'^ IJ>~ samples separately, obtaining 
B{B- ^ X{-i872)K-) X B{X{3872) -> J/V'tt+tt-) = 
(1.94±0.62) X 10-5 and (0.52±0.46) x 10"^ respectively, 
consistent within 1.8 standard deviations. 

The decay of a charmonium state into pJ/ip is a 
strongly suppressed isospin- violating process. In order to 
investigate the nature of the X(3872) state, we plot the 
invariant mass of the 7r+ tt^ system in both the X(3872) 
and the '0(25') region (Fig. ^. In the i){2S) case, the 
events are concentrated near the kinematic limit. Such 
behavior is not excluded for the X(3872), but the statis- 
tics are too small to allow a clear conclusion. Measuring 
both the TO^+^-and angular distributions with signifi- 
cantly greater statistics would provide important infor- 
mation on the nature of the X(3872). 

The search for the he is performed with the same 
frequentist technique in a ±10MeV/c^ mass window 



centered on m 

N, 



j/ip-R 






3526MeV/c2 |g. With Ndata = 9, 
comb = 6.9 ± 3.5, Npeak = 0.6 ± 1.5, and assuming the 
same efficiency e^, = (92 ± 1)%, we set a 90% C.L. limit 
B{B- -^ heK-) X B{he -^ J/ipTr+ir-) < 3.4 x lO'^. 
The probability that we would see a signal as large as 
the one observed from background fluctuations alone is 
39%. 

Finally, we search for B^ -^ J/tpD'^TT^ decays with 
D° -^ K-TT+. The decay D^ -^ K-tt+ would have 
an r.m.s. width of 5.4MeV/c^ in mx-TT+ as de- 
termined from Monte Carlo. We study this distri- 




FIG. 4: Distribution of mi^-^+ in events 

B^ -^ J/ip K~ ■k'^-k~ , with combinatorial background re- 
moved. Overlaid is an exponential fit. The arrow indicates 
the 3(T region expected for Z)° — > K~ir^ . 



bution in the same way we studied ttij/^^t^. The 
mx~Tr+ combinatorial-subtracted distribution (Fig. 0)) 
shows no significant structure, and it is therefore used 
to set a limit. We fit the background from other 
B~ — > J/^ K~TT^TT~ decays with an exponential func- 
tion of rnK-T!+ and obtain Np^ak = 2.9 ± 1.4. The 
frequentist approach described above, with Ndata = 10, 
Ncomb = 7.8 ± 2.8 and e/e^(2S) = 1.00 ± 0.07 yields the 
90% C.L. limit B{B- -^ J/iPD°t:-) < 5.2 x 10"^ This 
upper limit rules out the explanation of the inclusive 
J/ip momentum spectrum with intrinsic charm proposed 
in 0. 

In summary, we measured B{B^ —f J/ip K^ir^ir^) ~ 
(116± 7(stat.) ±9(sysi.)) X 10~^ with an error almost a 
factor two smaller than the present average iljj and we 
confirmed the observation of B- -^ X{3872)K- [IQ. 
We measured B{B- -* X{3872)K-) x S(X(3872) ^ 
J/V-TT+TT-) = (1.28 ± 0.41) X 10-5 and mx(3872) = 
3873.4 ± 1.4MeV/c2. We also studied the mj/^^„ and 
iT^K--K+ distributions searching for B^ —^ hcK~ and 
B~ — + J/ipD'^TT" decays and set limits on their branch- 
ing fractions, B{B^ -^ hcK^) x B{hc -^ J/i/jTr+Tr^) < 
3.4 X 10-^ and B{B- -^ J/^jD°ti-) < 5.2 x 10"^ at 90% 
C.L. 
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